Effect of linear density of states on the quasi-particle dynamics and small 
electron-phonon coupling in graphite 
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We obtained the spectral function of very high quaUty natural graphite single crystals using angle 
resolved photoelectron spectroscopy (ARPES). A clear separation of non-bonding and bonding bands 
and asymmetric lineshape are observed. The asymmetric lineshapes are well accounted for by the 
finite photoelectron escape depth and the band structure. The extracted width of the spectral 
function (inverse of the photohole life time) near the K point is, beyond the maximum phonon 
energy, approximately proportional to the energy as expected from the linear density of states near 
the Fermi energy. The upper bound for the electron-phonon coupling constant is about 0.2, a much 
smaller value than the previously reported one. 

PACS numbers: 74.25.Jb, 63.20.Ls, 79. 60.-1 
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Recent discoveries of novel physical properties in 
carbon-based materials such as superconductivitj^i^i^'^ 
and massless Dirac Fermions^ brought renewed interest 
in the electronic structure of graphite^'^. The peculiar- 
ity of the electronic structure of graphite has two as- 
pects: graphite is extremely two dimensional and is a 
semi-metal. These facts make it fundamentally interest- 
ing to study how dimensionality affects the dynamics of 
the doped carriers and how the carriers in graphite in- 
tercalated compounds (GICs) couple to the mediating 
bosons. In fact, there is a long standing issue in regards 
to the carrier dynamics in graphite, that is, whether 
the carriers are Fermi- liquid-like or not. This question 
motivates experimental studies of electronic structures 
of these materials by using, for example, angle-resolved 
photoemission spectroscopy (ARPES) and one can find 
a long history in the ARPES studies on graphite sin- 
gle crystalsi^i^iiSiiiii^ii^iiiii^. In addition, studies of 
graphite-related materials such as singlei^ and bilayer 
graphen&ii and GlC o-'^^d^ can be found. 

High quality ARPES data from graphite is difficult 
to obtain despite graphite's two dimensional, inert na- 
ture. The problems associated with ARPES experiments 
on graphite are due mostly to difficulty in proper sur- 
face preparation and to some extent to the low quality 
of the single crystals. For example, the extreme two- 
dimensional nature of graphite inevitably produces small 
flakes (many are small enough to be seen only under mi- 
croscopes) on the cleaved surfaces which ruins the mo- 
mentum resolution in ARPES. Such difficulties prevented 
one from obtaining good quality data to extract reli- 
able information on the many-body interactions such as 
electron-phonon coupling (EPC). Therefore, unless such 
difficulties are overcome, reliable information on many- 
body interactions can not be extracted from the data. 
As a result, the experimental data in regards to the elec- 
tron lifetime have been obtained mostly by time-resolved 



photoelectron spectroscopy on highly oriented pyrolytic 
graphite22i2i. 

Motivated by the renewed interest in the carrier dy- 
namics in graphite, we have performed ARPES studies 
on graphite single crystals. Our goal was to extract reli- 
able quantitative information on the EPC from ARPES 
data. To overcome the above mentioned difhculties, we 
exploited micro-spot ARPES and high quality natural 
graphite single crystals, and were able to obtain ARPES 
data with high enough quality to reliably extract quan- 
titative information on the electron-phonon coupling. In 
addition, we successfully applied the finite photoelectron 
escape depth and band structure effects in graphite to 
the analysis of the ARPES spectral functions for the first 
time (to our knowledge) . This allowed us to obtain an ac- 
curate measure of the EPC in graphite by using ARPES. 

ARPES experiments were performed at the beam line 
3A1 of the Pohang Accelerator Laboratory (preliminary) 
and 7.0.1 of the Advanced Light Source (main results). 
High quality natural graphite single crystals with sizes 
larger than 1 cm were cleaved ex situ. The samples were 
subsequently introduced to the vacuum chamber and an- 
nealed for half an hour at about 900° C in a vacuum better 
than 6.0x10"^° Torr. ARPES measurements were per- 
formed at 25 K with an energy resolution of about 40 
meV which mostly stemmed from 0.0075 momen- 
tum resolution in combination with the fast dispersion 
of the bands. The pressure was better than 5.0x10"^^ 
Torr. For the experiment, we exploited the small beam 
spot (« 40//m) to locate a flat region without flakes. A 
photon energy of 85 eV was used so that was near K. 
Electronic structure calculation was done by using the 
SIESTA code^^ based on pseudo-potential method. 

Fig. 1(a) shows the calculated electronic band disper- 
sion of graphite along the high symmetry lines, F-M-K-F 
(solid lines) and A-L-H-A (dashed). Also shown in the 
figure as an intensity map is the experimental data with 
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FIG. 1: (a) Comparison of experimental data and calculated 
band structure of graphite along the high symmetry lines, 
(b) Intensity map at a constant binding energy of 1.75 eV. 
Brillouin zones are drawn as solid lines. 



the photon energy tuned for kz =K. Other than the fact 
that the calculated band structure has to be expanded, 
which is well knowo^.^ theoretical and experimental re- 
sults match well. The two bands near the Fermi energy 
are Pz-derived tt bands and are split due to the inter-layer 
interaction, resulting in anti-bonding (AB), non-bonding 
(NB) and bonding (BB) bands. Among them, NB and 
BB are occupied and can be seen by ARPES. The split- 
ting UJB between the NB and BB is about 0.7±0.1 eV^. 
The data can also be plotted in momentum space at a 
constant binding energy (1.75 eV) as shown in Fig. 1(b). 
The plot shows an apparent hexagonal BouUiouin zone 
expected for graphite. It is also clearly seen that rounded 
triangles around the K points are split. The bands have 
approximate linear dispersion with point-like Fermi sur- 
faces at the K points, resulting in an approximate linear 
electronic density of states near the Fermi energy (linear 
DOS). 

Fig. 2(a) shows high resolution data taken along the 
M-K-F line as marked in Fig. 1(b). The two peaks are 
sharp and clearly split. We note that there are no defect- 
induced states as reported in Refj^. In addition, the 
background at high binding energies is very small, in 
strong contrast to the previously reported dat a^'^'^°i^^'^^ . 
We attribute these to the high quality of the crystals used 
in our experiment. Fig. 2(b) shows the energy distribu- 
tion curve (EDO) at the K-point (arrow in panel (a)). 
One peculiar aspect of the data is that the region be- 
tween the two peaks are somewhat filled up. As a result, 
the BB has a tail on the lower binding energy side con- 
trary to the usual case. This aspect of the data can be 
understood as follows. Even though we tuned the photon 
energy to the K point, the finite escape depth of the pho- 
toelectron introduces an uncertainty in kz, Ak^ = l//i 
where fj, is the escape depth. This uncertainty in de- 
picted in Fig. 2(c) necessarily brings in broadening due 
to the kz dispersion, which is known as the photoelec- 
tron life time. Note that the NB (BB) has its dispersion 
maximum (minimum) at fc^ = and as a consequence 
the broadening effect will be one-sided towards the band 
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FIG. 2: (a) High resolution ARPES data taken along the M- 
K-F as marked in Fig. 1(b). (b) A typical EDC (circles), taken 
at the k value marked in panel (a) and the fit (line) using the 
model discussed in the text, (c) Schematic of the disper- 
sions along the kz direction for non-bonding (upper curve) 
and bonding (lower curve) bands. Due to the finite escape 
depth {/J,), a range of kz contributes to the spectral functions. 
The contribution probability of each k^ is represented by the 
height of the graph, (d) Resulting spectral function near K 
point without lifetime broadening. 



position at kz — n/c where the two bands collapse. Fig. 
2(d) depicts a model spectral function when all these ef- 
fects are accounted for. Only when such effects are con- 
sidered can one reliably extract the lifetime broadening. 
We used fj, = 7 A for the fitting^^ and the model function 
in Fig. 2(d) is convolved with a Voigt function with the 
Gaussian width set to the total energy resolution of 40 
meV. A typical fit is shown in Fig. 2(b) as the solid line. 

Before we look into the low energy dynamics of 
graphite, it is necessary to discuss how the couplings to 
phonons and other low energy excitations are affected by 
the linear DOS of graphite. The EPC in metallic systems 
is well understood^® and has been extensively studied re- 
cently by using ARPES^li^^. However, there is a relative 
lack of understanding of EPCs in semi-metals. In fact, to 
our best knowledge, EPC effects on the ARPES lineshape 
for semi-metals with linear DOS have not been discussed. 
As a consequence, models derived for metals have been 
improperly used in graphite studies^. 

To understand the EPC in a system with a linear DOS, 
we evaluate the Feynman diagram for a single phonon 
process shown in Fig. 3(a). We first look at a single 
band system for simplicity and will extend the analysis 
to the graphite case. The contribution from the process 
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to the imaginary part of the self energy S"p due to EPC 
at r=0 if^^ 

Kp{uJk) = XI / l< k',q\Hi\k,0 >\^ f{-ujk')S{uJk-(^k'-uJ„, 

(1) 

where k and fc' are crystal momenta of the holes and / 
Fermi-Dirac function at T = 0, respectively. The phonon 
momentum q — k — k\ and v represents the phonon 
mode. tl\ is the Hamiltonain for the EPC, and it can be 
written as 



(2) 



where g is the transition probability amplitude for partic- 
ular k and q. Assuming Einstein phonons and constant 
9, 



< k',q\gcl_^Ckbl\k,0 >\^S{uJk - ^k' - U!„^q)dk' 



/ 5{ujk - Wfe' - U!u,q)dk' = g'^T>{ujk ~ '^v) 



(3) 



where T) is the electronic DOS. The expression for S^'p 
then becomes 



(4) 



For a metal, 2? is approximately constant in the imme- 
diate vicinity of Ep and above equation reduces to the 
expected step function. On the other hand, for a lin- 
ear DOS, S"p(ti;) increases linearly beyond the phonon 
energy as shown in Fig. 3(b}^. 

We may now extend the analysis to the two band case. 
The EPC Hamiltonian becomes 



^1 = ^9iA- 



1,3 



Ck,i 



(5) 



where i and j are band indices with i representing the 
initial band, gij with i = j is the EPC for intra-band 
transition while gij with i ^ j for inter-band transition. 
The result is shown schematically in Fig. 3(c)"^°. Note 
that E"p for the NB has an upward kink at the energy 
where the BB starts. The kink and slopes of the self 
energy curves carry the information on the EPC matrix 
element gij. 

The difficulty in extracting the EPC constant A comes 
from the fact that it is almost impossible to guess the 
bare bands. In metals, the bare band can be approxi- 
mated by a linear band (possibly with an additional small 
parabolic term) over the energy window of interest which 
can be obtained by fitting the experimental dispersion. 
This method, however, cannot be used for semi-metals 
and insulators for an obvious reason. In fact, the NB 
band of graphite near K does not form a Dirac cone 
but has a fairly strong curvature as is well known from 
the band calculation. Therefore, the standard method of 
obtaining the real part of the self energy E' by taking 




q=0 mode 



FIG. 3: (a) Feynman diagram for the EPC under consider- 
ation, (b) Schematic of the E" for the phonon decay for a 
single band system with a linear DOS. (c)E" for a double 
band system, (d) Other possible decay channels for photo- 
holes, red for decay through electron-hole creation and blue 
for g=0 mode decay. Note that k' and kh can be changed for 
the electron-hole creation case. 



the difference between the bare band and the experimen- 
tal dispersion can not be used. Alternatively, E' can be 
obtained from E" through the Kramers-Kronig relation. 
Then the EPC constant A can be extracted from the real 
part of the self energy using the following formulap2£. 



(6) 



We point out that this method usually requires higher 
quality of data compared to the standard method to ex- 
tract the similar quality of information. 

Finally, we discuss other decay channels (see Fig. 3(d)). 
The photo-hole may decay by creating an electron-hole 
pair or by exciting a q—Q mode in addition to emitting 
a phonon. For the decay through electron-hole pair cre- 
ation, the energy and momentum conservations mandate 
Aa;=a;fe' — LJk=u!e — tUh and Ak=k -k^ke-k/j. Note that 
both k' and k are in the same (lower) Dirac cone while kg 
and k/i are in different Dirac cones. In such case, for the 
equations to be met, both k'-k and ke-kh vectors should 
lie in the steepest decent line on the 2 dimensional coni- 
cal dispersion (the configuration illustrated in the figure) 
and the available phase space volume for the transition 
in Acj-Ak space is zero^. For the q = mode decay, the 
contribution is considered to be photo-hole momentum- 
independent as the two bands are almost parallel to each 
other. Note that NB band is not susceptible to g = 
mode decay as one can see in the figure, giving a con- 
stant upward shift only to E" of the BB. 

With the finite escape depth considered, the half 
width of the Lorentzian is extracted and plotted in Fig. 
4(a). Filled and empty symbols represent the BB and 
NB, respectively as marked in Fig. 2(a). Overall, the 
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FIG. 4: (a) E" (half of the Lorentzian width) vs. binding 
energy, (b) Calculated PDOS of the NB and BB. (c) E" for 
the NB over the entire energy range range by combining the 
experimental data and the calculated PDOS. (d) E' obtained 
by Hilbert transforming the E" in panel (c). 

width monotonically increases without particular fea- 
tures. There are smooth slope changes at around 0.4 
eV which roughly coincide with energy where kz disper- 
sion of the NB saturates. Therefore, we attribute this 
to the yet unaccounted band structure (in combination 
with the finite escape depth) effect"^^. When such an 
effect is considered through additional modeling, we ex- 
pect the width will increase more linearly. While the 
values for the NB, when converted to life time, are rea- 
sonably compatible with the values obtained by time re- 
solved techniques^i^l, they are much smaller than the 
previously reported values^. 

From the measured widths, we find the the follow- 
ing additional facts. Firstly, the fact that the observed 
widths do not show a high order behavior (~ w^) re- 
veals that the EPC is the dominant decay channel for 
the photo-holes. This is expected because the available 
phase space is very small as discussed above. Secondly, 
no (or very weak at most) kinky feature at the optical 



phonon energy 0.2 eV) is observed. It was previously 
reported that a strong EPC was observed in the imag- 
inary part^ of the self energy. Our observation shows 
it is not the case (at least near the K point). Thirdly, 
the inter-band decay, if exists, is much weaker than the 
intra-band evidenced by the lack of an upward kink in 
the NB band curve at the onset energy of BB (see Fig. 
3(c)). This allows us to have a convenience of treating the 
NB and BB independently. As will be clear, this makes 
the task of calculating the EPC constant much simpler. 
Fourthly, as the inter-band decay is very small, the width 
of the BB at the onset energy can not be explained by 
the EPC. It has to come from another mechanism such 
as the q = mode decay as discussed earlier which is 
not well-understood. Lastly, the slopes of the curves are 
somewhat fc-dependent (that is, the slopes for the K-M 
cut are steeper). This fc-dependence suggests that there 
is fc-dependent EPC in the system. Such behavior is in- 
deed predicted theoretically — and is also observed in the 
study of a single layer graphene^. 

Even though the constant offset in the E" of the BB is 
not well-understood, one can estimate the EPC of the NB 
as the inter-band transition is weak. As discussed above, 
calculating the EPC constant A requires knowledge of S" 
over the entire energy range. As it is not available, we 
approximated S" at high energies by the scaled, calcu- 
lated partial DOS (PDOS) to match the the experimen- 
tal data, noting that S" is proportional to PDOS as in 
equation (3). The resulting E" for the NB along the K- 
M direction is plotted in panel (c). E' is then obtained 
by a Hilbert transformation of the symmetrized E" and 
plotted in panel (d). The Xnb value obtained from E' 
is about 0.12, which is much smaller than the previously 
estimated value from ARPES^. For the K-F direction, 
A was 0.2. Therefore, we set 0.2 as the upper bound for 
EPC constant A. This value is very small compared to A 
values in typical metals as expected for semi-metals with 
very small DOS at the Fermi energy. 
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This file contains supplementary information on the 
the self energies due to the electron-phonon coupling 
(EPC) and electron hole pair creation in graphene and 
graphite which have characteristic linear electronic den- 
sity of states near the Fermi energy (linear DOS). We first 
investigate the scattering process through phonons by 
using schematic diagrams for the single (graphene) and 
double (graphite) conical band structures. It allows one 
to intuitively understand how characteristic linear DOS 
contributes to the scattering rate of photo-hole. Then 
we discuss the scattering process by electron-hole pair 
creation by considering the available phase space for the 
process. We show that this process gives little contribu- 
tion to the scattering rate of a photo-hole. 



I. SCATTERING RATE DUE TO 
ELECTRON-PHONON COUPLING 

A. Single band case 

In this case, we assume a conical band structure with 
the Fermi energy at the apex of the cone. There also is an 
unoccupied conical band above the Fermi energy, forming 
a Dirac-cone-like band structure as shown in Fig. 1(a). 
Assume an Einstein phonon at energy loq. If a photo-hole 
with momentum k (filled circle) is created by a photon 
as shown in Fig. 1(a), it can be filled by an electron 
with energy of ujk' = — and momentum k' = k — q 
(empty circles) where q is the phonon momentum. We 
further assumed for simplicity the scattering amplitude 
g is momentum independent. In that case, the scattering 
rate is proportional to the number of such k' states, thus 
the DOS at ojk ~ ojq- Note that if the binding energy of 
the photo- hole is smaller than the phonon energy luq, the 
scattering cannot occur because there are no electrons 
with sufficient energy to emit a phonon. Therefore, the 
self energy of photo-hole as a function of the binding 
energy is proportional to V^uJk — ujo) and looks like the 
schematic shown in Fig. 1(b). 



B. Double band case 

For the double band case, another conical band exists 
ujb below the Fermi energy as shown in Fig. 2(a). One 
needs to investigate the cases where the initial photo- 
hole is in non-bonding and bonding bands separately. 
We again count the the number of available electrons 
that can fill the photo-hole. We first consider the case 
in which the photo-hole is created in non-bonding band. 
If the binding energy of the photo-hole is smaller than 
ujb+'-^Oj the decay rate (the number of electrons that can 
fill the photo-hole) will be the same as that in the single 
band case. When the binding energy becomes larger than 
(jJb+ojq, the additional electrons (triangles in Fig. 2(a)) 
in the second band begin to contribute to the scatter- 
ing process. This additional contribution is proportional 
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FIG. 1: (a) Schematic diagram of the photo-hole decay by 
emitting a phonon in the single band case, loq is phonon 
energy and k and k' are initial and final photo-hole states, 
respectively, (b) Schematic of the resulting decay rate (the 
imaginary part of self energy) vs binding energy. It is zero up 
to ujQ and then increases linearly. 
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FIG. 2: (a) Schematic diagram of the photo-hole decay in the 
double band case, ujb is the energy of the apex of the lower 
conical band, (b) The imaginary parts of the self energies vs 
binding energy for photo-holes created in non-bonding and 
bonding bands. 



to the partial DOS of the second band at Uk — i-^b and 
thus linearly increases. Note that the additional slope 
of the linear increase depends on the inter-band scatter- 
ing amplitude gtj (i ^ j). The resulting decay rate is 
schematically shown in Fig. 2(b). We now investigate 
the case for a photo-hole in the bonding band. If the en- 
ergy of the photo-hole is larger than lob but smaller than 
lub +^^0; the photo- hole can be filled with only electrons 
in non-bonding band via inter-band transition. Mean- 
while the photo-hole can decay via inter- and intra-band 
transitions when the energy is larger than ujb + loq. The 
resulting decay rate is also schematically shown in Fig. 
2(b). If the electron-phonon couplings are different for 
bonding and non-bonding bands, the overall slopes of 
self energies will be different as shown in Fig 2.(b). 
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FIG. 3: (a) Schematic of electron-hole creation(ke — k^) and 
transition of a photo-hole (ke — k') in a conical band structure, 
(b) Available phase spaces for k^ — kh (hatched region) as well 
as kg — k' (grey region) processes. Two regions overlap only 
at the boundary. 



a transition to a k' state by creating an electron-hole 
pair (denoted as ke and kh, respectively in Fig 3. (a)). If 
we plot the energy difference {Auj=ujk — cok') as a func- 
tion of momentum difference Ak=k — k', the possible 
transitions occupy the lower right of the Aw-Ak phase 
space as shown in Fig. 3(b) (grey). In a similar way, 
we can go through the same process for the electron- 
hole pair creation. One can show that Au}=u!e — ujh and 
Ak=ke — kh of the electron-hole creation processes oc- 
cupy the upper left of the phase space (hatched). We 
note that only the boundary is shared by two regions. 
Energy and momentum conservations enforce the decay 
through the electron-hole pair creation to occur in the 
phase space points that belong to both regions. There- 
fore, the available phase space volume for the decay is 
zero and thus the transition rate. 



II. AVAILABLE PHASE SPACE FOR 
ELECTRON-HOLE PAIR CREATION 



A photo-hole with the momentum k is created in the 
Dirac cone as shown in Fig 3. (a). It may go through 



